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Abstract 

A class of Poisson embeddings of reduced, finite dimensional symplectic vector spaces 
into the dual space ~q* r of a loop algebra, with Lie Poisson structure determined by the 
classical split i?-matrix R = P+ — P_ is introduced. These may be viewed as equivariant 
moment maps inducing natural Hamiltonian actions of the "dual" group <3r = © + x <5~ 
of a loop group on the symplectic space. The i?-matrix version of the Adler-Kostant- 
Symes theorem is used to induce commuting flows determined by isospectral equations of 
Lax type. The compatibility conditions determine finite dimensional classes of solutions 
to integrable systems of PDE's, which can be integrated using the standard Liouville- 
Arnold approach. This involves an appropriately chosen "spectral Darboux" (canonical) 
coordinate system in which there is a complete separation of variables. As an example, 
the method is applied to the determination of finite dimensional quasi-periodic solutions 
of the sine-Gordon equation. 



0. Introduction. 

In [AHP, AHH1-AHH3] a unified framework was developed, describing a wide class 
of integrable finite dimensional Hamiltonian systems, as well as finite dimensional solutions 
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ported in part by the Natural Sciences and Engineering Research Council of Canada and the Fonds FCAR du 
Quebec. 
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to integrable systems of PDE's. The basic idea is to represent all these systems in terms of 
commuting isospectral flows determined by Lax equations within finite dimensional Poisson 
subspaces of loop algebras consisting of orbits that are rational in the loop parameter. The 
commuting flows are generated by Hamiltonians that are spectral invariants, or equivalently 
A<i*-invariants on the space q + * dual to the half of the loop algebra admitting holomorphic 
extensions to the interior disk of a circle within the complex plane of the loop parameter. 
According to the Adler-Kostant-Symes (AKS) theorem, such invariants generate commuting 
flows with respect to the Lie Poisson structure on g+* and are determined by equations of 
Lax type. 

Finite dimensional, generically quasi-periodic solutions to integrable systems of PDE's 
arise within this framework as the compatibility conditions satisfied by the matrix elements, 
within a given representation. The procedure developed in [AHP, AHH1-AHH3] consists 
of three steps; first, a suitably defined moment mapping is used to embed symplectic vector 
spaces quotiented by Hamiltonian group actions as Poisson subspaces of the space The 
image consists of rational coadjoint orbits with fixed poles. The AKS theorem is then used to 
define the Hamiltonians and determine the dynamical equations, with the ring of invariants 
generated by the invariant spectral curve given by the characteristic equation of the loop al- 
gebra element. The generic systems so obtained may be subjected to further reductions under 
both discrete and continuous Hamiltonian symmetry groups. Finally, a suitably defined Dar- 
boux (canonical) coordinate system is introduced, associated to the divisor of zeroes of the 
matrix representating the loop algebra element on the spectral curve. Within this "spectral 
Darboux" coordinate system, there is a complete separation of variables which, through ap- 
plication of the Liouville-Arnold integration method on the invariant Lagrangian manifolds, 
reduces the problem of integration to quadratures given by Abelian integrals. The resulting 
flows are thus seen to be linearized by the Abel map in the natural linear coordinate system 
on the Jacobian variety of the spectral curve, or some quotient thereof. Through the Jacobi 
inversion method, the matrix elements may be expressed in terms of theta functions, thereby 
producing, through classical Hamiltonian methods, the types of solutions typically obtained 
through more sophisticated algebro-geometric means [AvM, KN, D, AHH1] . 

This approach has been applied to a wide variety of integrable systems, both finite and 
infinite dimensional, such as: the cubically nonlinear Schrodinger equation and various multi- 
component generalizations thereof [AHP, AHH3, HW1, Wl], the sine-Gordon equation 
[HW2, TW], the massive Thirring model [W2] the Neumann and Rosochatius systems 
[AHP, AHH1, AHH3] and various generalized tops [AHH2, HH]. 

In the subsequent section, an extended form of the moment map embedding method of 
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[AHP, AHH2] will be given, based on the classical .R-matrix technique, in which the two 
halves of the loop algebra g + : g~ are treated on an equal footing. In section 3, the process of 
discrete reduction, together with the Liouville- Arnold linearization method, will be illustrated 
for the case of the sine-Gordon equation. Full details for the latter may be found in [HW2]. 



1. Moment Maps to Loop Algebras and the Classical R— Matrix. 

We first establish notations for loop groups and algebras. Let (M, uj) denote the sym- 
plectic vector space whose elements are pairs (F, G) of N x r matrices 

M = {(F, G) e M N,r x M N,r }, (1.1) 

with symplectic form: 

u = tr (dF T AdG). (1.2) 

Let (25 denote the group of smooth loops in 0l(r) (real or complex), or some subgroup thereof, 
viewed as invertible matrix- valued smooth functions g(X) defined on a smooth, simple closed 
curve T enclosing the origin of the complex A plane. Let <S + , <3~ be the subgroups of loops 
admitting holomorphic extensions, respectively, to the interior and exterior regions r + , T~ 
(including oo), such that for g e (J5 _ , g(oo) = I. The Lie algebras corresponding to 35, (S^, 
are denoted g and ^ respectively. Their elements are smooth gt(r) valued functions on T, 
with elements X + e Q + admitting holomorphic extensions to r + and X_ E g~ to T - , the 
latter satisfying X(oo) = 0. 

We identify g as a dense subspace of the dual space q* through the pairing: 

<X,y>:=^/ r tr OT (A))f, 
X eg* , Y eg. (1.3) 

Using the vector space decomposition 

= 0"©0 + , (1-4) 
this allows us to identify the dual spaces as 

0-*~0+, (1.5) 

where g_, denote, respectively, the subspaces of g consisting of elements admitting a 
holomorphic extension to T~ or T + , with X + e g + satisfying X + (0) = 0. 
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Let 

p± •■ — > 

1.6 

P ± :Ii — ► X ± 

be the projections to the subspaces g^ 1 relative to the decomposition (1.4) and define the 
endomorphism R : g — > g as the difference: 

R:=P + -P_. (1.7) 

Then is a classical R-matrix [S], in the sense that the bracket [ , ]_r defined on g by: 

[X,Y]*:=^X,Y] + i[X,i?Y] (1.8) 

is skew symmetric and satisfies the Jacobi identity, determining a new Lie algebra structure 
on the same space as g, which we denote g^. The corresponding group is the "dual group" 
(5r = <5~ x (J5 + associated to (5. The adjoint and coadjoint actions of (&r on g R ~ g* R are 
given by: 

Ad R (g) : (X+ + XJ) .— g+X+gl 1 + g-X.gZ 1 (1.9a) 
Ad^) : (Y+ + y_) .— > (g.X +9 Z 1 )+ + (g+X.gl 1 )., (1.9b) 

^±G0 ± , Y ± eg±, g±e® ± . (1.9) 
The Lie Poisson bracket on g^ ~ g* dual to the Lie bracket [ , ]r is: 

{f,g}\x=< [df,dg] R ,X> (1.10) 

for smooth functions /, g on g^ (with the usual identifications, df\x, dg\x G g* ~ g). 

Let A G gl(iV) be a fixed N x N matrix having no eigenvalues on T, and define the 
following action of <3r on M 

& R :M — > M 
g(X):(F,G)^(F g ,G g ) (1.11) 

where (F g , G g ) are defined by: 

Fg := F ~ i £ (A ~ A/) ~ lF ^+ 1(A) ~ ^ zl(A)) dX (L12a) 

G g :=G - £(A T - Xiy'G (gl(X) - £(A)) dA. (1.12b) 
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It is straightforward to verify, using the identity 

A — a 

and residue calculus, that the (&r composition rule is satisfied, so (1.11), (1.12a,b) does, 
indeed define a 0^-action on M. A similar calculation shows that this action preserves the 
symplectic form (1.2) and is, in fact, generated as a Hamiltonian action by the equivariant 
moment map: 

J a ,y :M — > q* r 
J a ,y(F, G) = XY + XG t (A - XI)~ 1 F, (1.14) 

where Y e 0l( r ) is any r x r matrix. The constant term Ay may be included without affecting 
the equivariance of the moment map (1.14) since Ay is an infinitesimal character for the Lie 
algebra g R : 

<XY,[X,Y} R >=0, VX,Yeg R . (1.15) 

The splitting Ja,y(F, G) = J+(F,G) + J-(F,G), J±(F,G) G g± is determined by the pole 
structure, with the Ay term plus the poles at eigenvalues of A in T~ included in the ~ 0+ 

part and the poles at eigenvalues in r + in the (0 + )* ~ 0- part. 

The fibres of the Poisson map Ja,y are the orbits of the stability subgroup Ga '■= 
Stab(A) C Gl(N) under the Hamiltonian Gl(N) action defined by 

g : M — > M 

g:(F,G)^(gF,(g T )- 1 G). (1.16) 

On a suitably defined open, dense set, this action is free, allowing us to define the quotient 
Poisson space M/Ga ■ Since the map Ja,y '■ M — > g R passes to the quotient, defining a 
1-1 Poisson map Ja,y '■ M/Ga — > 0^, we may identify M/Ga with the image Poisson space 
Qa c &*r consisting of elements Af(X) of the form 

n m 
i=l i i= i ^ A UL% > 

where {ai}i=i j ... n are the eigenvalues of A, {ii}i=i,... n are the dimensions of the corresponding 
Jordan blocks of generalized eigenspaces, and the ranks and Jordan structures of the r x r 
matrices iVi^ are determined by the multiplicities of the eigenvalues of A and its Jordan 
structure. 
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The Hamiltonian equations on the phase space g\ ~ Mj Ga are generated by elements of 
the ring Z\ consisting of Ad* -invariant polynomials on g*, restricted to the Poisson subspace 
C g|j. According to the Adler-Kostant-Symes (AKS) theorem, in its .R-matrix form [S], 
the elements of this ring Poisson commute, generating commuting Hamiltonian flows, and 
Hamilton's equations for e X\ have the Lax form: 

^^ = [(#)+,a/] = -mu at], 

(d(f>)±e9 ± , #=(#)+ + (#)- (1.18) 
This implies that the spectral curve <S, with affine part defined by the characteristic equation 

det(£(A) -zl) :=V(X,z) =0, (1.19) 

where 

f.(\) -= a - x 



£( A ) : = ^AT(A), a(A) := f[(X - a*)"* (1.20) 



i=i 

Y 



is invariant under these flows, and its coefficients generate the ring X A . 

Integrable systems of PDE's then arise as the equations satisfied by the matrix elements 
of M implied by the compatibility conditions 



d(d(p) 



±-<^l± + [d(0) + ,d(^ + ]=O, (1.21) 



dx dt 

where t, x are the respective flow parameters for the Hamitonians and -0. The flows may 
be integrated through a standard algebro-geometric construction that leads to a linearizing 
map to the Jacobi variety of the spectral curve S ([KN, D, AHH1]). To obtain interesting 
examples, the generic systems so obtained must usually be further reduced by certain contin- 
uous or discrete symmetry groups. The easiest way to arrive at the linearization involves the 
introduction of a suitably defined "spectral Darboux" (canonical) coordinate system associ- 
ated to the divisor of zeroes of the eigenvectors of C(X) on the spectral curve. The general 
construction is detailed in [AHH3] . Its application to the particular problem of determining 
finite dimensional quasi-periodic solutions to the sine-Gordon equation is described in the next 
section. Full details for this case may be found in [HW2]. 
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2. Quasiperiodic Solutions of the Sine-Gordon Equation 

To obtain the sine- Gordon equation 

d 2 u d 2 u 



dx 2 dt 2 



= sin(u) (2.1) 



we shall need commuting isospectral flows in the twisted loop algebra iu(2) := 5U ( ' 1 ' ) (2). This 

is the subalgebra of 0((2) consisting of elements X(X) that are invariant under the three 
involutive endomorphisms: 

(7i : X(X) i — ► JX T (X)J (2.2a) 

a 2 :X(X)i — >-X^(X) (2.2b) 

£7 3 : X(X) i — ► tX(-X)t (2.2c) 

where 

r:=(i. °X "iV (2-3) 



v o — i y ' V 1 °, 

The first of these implies that X is traceless, the second that it is in su(2) and the third that 
it is in the "twisted" subalgebra su(2). 

We choose the matrix Y as 

Y = J=(° 1 "J). (2.4) 

T as a circle centred at the origin of the A-plane, and A as a diagonal 2N x 2N dimensional 
matrix, having distinct eigenvalues, all inside T + , of the form 

A = diag (oji, ai,..., a p , a p , -a*!, . . . , -a p , ifop+i, -ifap+i, • • • , 0n, -0n), (2.5) 

where {aj}j=i,..., p have nonvanishing real and imaginary parts and {(3j}j=2p+i,...,N are real. 

The space M = {(F,G) G M 2N > 2 x M 2N > 2 } must correspondingly be reduced by the 
endomorphisms 



E 



E 2 
S 3 



(F, G) i — ► (GJ, -FJ) (2.6a) 
^ (-JG,JF) (2.6b) 
(F, G) i— > (-ilCFr, iKGt) (2.6c) 



where 



J" = diag (J, J,..., J) 



(2.7) 
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is the block diagonal 2N x 2N matrix with 2x2 blocks equal to J and 



/ 



<-2p 



-I 



2p 



J 



\ 



\ 



J/ 



(2.8) 



is the 2N x 2N matrix consisting of a Ap x Ap block formed from I 2p , the 2p x 2p identity 
matrix and its negative, and N — 2p diagonal 2x2 blocks J. 

The moment map J\ of eq. (1.14) then intertwines the automorphism groups generated 



by £i, £ 2 , £3 and by <7i, cr 2 , 03. Denoting by 



F 2 j-i 



'2j 



G2j-1 



the consecutive pairs 



of 2 x 2 blocks in (F, G), the fixed point set Ms C M under Ei, E 2 , E 3 consists of (F, G) with 
2x2 blocks of the form 



F 2j -i 

F2J 



F 



2j+2p-l 
7 2j+2p 



<2j 



<Pj 7j 

7j -Vj 
i(pj -%y 
Hj i<Pj 



G 2 j-i 
G 2j 

G 2 j+ 2p -i 
G 2 j+ 2p 

G 2j -i 
G 2j 



-y Vj 

<Pj 7j 

y m 
-iipj y 



3 = 

j = i, 



"7j *7j 
-*7j 7j 



We have 

E^c<; = u), E^cu = a), E3CU = u> 
so M s is a real symplectic space, with symplectic form: 

P N 

Q = A y^(c?7j A d<£j + <i7j A d^j) + 4z d^j A Gfy,-. 

j=l j=2p+l 

The restriction of J\ to M a , which we denote J, is given by 



J(F,G):=Af(\)=\r i Q - ) + 2X 



-1 



/6(A) -c(A)A 
^c(A) -6(A) J 



e su(2). 



where 



-Vili + fjlfj 



N 



i=i 



■w = E 

i=i 



a 2 -A 2 a 2 -A 2 



• 'A \ 

^ /3 2 + A 2 



N 



a? -A 2 a 2 - A 2 



i=2p+l •? 



N. 



(2.9a) 

(2.9b) 
(2.9c) 

(2.10) 
(2.11) 



(2.12) 



(2.13) 
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Now, choose Hamiltonians H^,H V e T\ as: 

Ht(X) = ~ jf tr (^(X(A)) 2 ) dA (2.14a) 

= ii tr ($^ (x(a))2 ) rfA ' (2 - i4b) 

where 

p JV 

a(A) = n[( A2 -«?)( A2 -«')] II ( A2 + ^)- (2- 15 ) 

j=l fc=2p+l 

Then the Lax form of Hamilton's equations may be written as 

dM 



dr\ 

where, setting 



->///,( A') .X (2.16a) 
r/A [^(A/) + ,Ay], (2.16b) 



C(\) = ^) A/(A) = A 2Ar_1 £ + A 2JV " 2 £i + • • • + £>2N—i + a(X)Y, (2.17) 
A 



we have 



dH^M)- = - T (C 2N -i + a(0)Y) (2.18a) 
A 

dH n (N) + = C + XY. (2.18b) 
The spectral curve is of the form 



det(£(A) - zl) = V(X, z) = z 2 + a(X)P(X) = (2.19a) 

P(A) = P + A 2 Pi + • • • + \ 2N - 2 P N -! + X 2N , (2.19b) 

with all coefficients Pi in and, in particular 

Ht = P , H r} = -P N _ 1 . (2.20) 

Choosing the invariant level set 

a(0)P o = det(£ 2JV _i + a(O)Y) = (2.21) 
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implies that 

£ 2JV _! + o(0)y = ^ _;_ m ^ J , (2.22) 
where u is real. Setting 

£ > = x + t, rj = x-t, (2.23) 

the compatibility equations for (2.16a,b) then imply that u satisfies the sine-Gordon equation 
(2.1). 

Since the spectral curve C determined by (2.19a,b) is invariant under the involution 
(z, A) h- > (z, —A), it is a two-sheeted covering of the hyperelliptic curve with genus N — 1 
whose affine part given is by 

z 2 + a(E)P(E) = 0, (2.24) 

where 

P(A 2 ):=P(A), a(A 2 ):=a(A), A 2 = E. (2.25) 
We also define functions 6, c by 

6(A 2 ) = 6(A), c(A 2 ) = c(A). (2.26) 
Setting z := zX, we obtain the genus iV hyperelliptic curve C with affine part given by 

i 2 + Ea(E)P(E) = 0, (2.27) 
which has two additional branch points at E = 0, oo. 

Following the general method of [AHH3], we define on C the divisor of degree N with 
coordinates (E^, C/j,)^=i,...,n determined by solving the equation 

2c(£ M ) + 1 = (2.28) 

for {£? M } At =i,...,jv and substituting in 



c »-« B^iy-T^- (2 - 29) 

This defines the zeroes of the eigenvectors of £(A) on the spectral curve C. In terms of these, 
we have: 



u 



-iln ^2 JJ £ M j + e7T (2.30) 



where e = 1,0 for AT even or odd, respectively. The {E^, C,^} ^=i,...,n may be viewed as 
coordinate functions on the coadjoint orbit through J\f = J e iu#(2), and the following result 
is key to the linearization of the flows (cf. [HW2] for details): 
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Proposition. The functions (E^, C^) h=i,...,n form a Darboux coordinate system on the coad- 
joint orbit passing through Af(X). The corresponding orbital symplectic form is 

N 

u = J2 dE n A d C^ = -dO. (2.31) 

This may be seen either by an explicit coordinate transformation from (2.11), by verifying 
the following expression for the canonical 1-form 

V N N 

9 = 4 j^iipid^i - 7id^) + 4i "fjdli = Yl ^ dE ^ ( 2 - 32 ) 

i=l j=2p+l M=l 

or as an application of the general "spectral Darboux coordinates" theorem of [AHH3] . 

It now follows directly from the definition of the coordinates (E^, Cn)fi=i,...,N that on 
the invariant level sets (Lagrangian manifolds) determined by {Pi = Cj}^ 1 , the 1-form 
@\{p i= ci} = dS may be integrated to yield the Liouville generating function 

s ^=££fM dK (2 - 33) 

Within the new canonical coordinate system {Qi, Pi}i=o,...N-i defined by 

«< = w = ~\ £ t i = ... iE ' < 2 - 34 > 



u=i JE ° \J-Ea(E)P(E) 

the flows for all the Hamiltonians in the spectral ring generated by the Pj's are then linear. 
In particular, integrating Hamilton's equations for H^, H v gives 

N r E » E i 

/ i = =dE = C i + 25 i ,ot-25 i , N - 1 ri. (2.35) 



„=i JB ° \J-Ea(E)P{E) 

This may be interpreted as the linearizing Abel map to the Jacobi variety of C. By a stan- 
dard inversion procedure (cf. [HW2], the function u may be explicitly expressed in terms of 
quotients of the associated theta functions: 

u = -2zm— — - — r — — + c, (2.36) 

Q{A{oo) - Urj - V£ - B - k) 

where the constant vectors A, U,V,E> e are obtained from those defined by the slopes 
and integration constants on the RHS of eq. (2.35) by applying the matrix that transforms 
the differentials appearing in the integrands on the LHS to a normalized canonical basis of 
abelian differentials. 
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